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SmO thin films: a flexible route to correlated flat bands with nontrivial topology 
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Using density functional theory based calculations, we show that the correlated mixed-valent 
compound SmO is a 3D strongly topological semi-metal as a result of a 4/-5(i band inversion at the 
X point. The [001] surface Bloch spectral density reveals two weakly interacting Dirac cones that are 
quasi-degenerate at the M-point and another single Dirac cone at the T-point. We also show that the 
topological non-triviality in SmO is very robust and prevails for a wide range of lattice parameters, 
making it an ideal candidate to investigate topological nontrivial correlated flat bands in thin-film 
form. Moreover, the electron filling is tunable by strain. In addition, we find conditions for which 
the inversion is of the 4/-6s type, making SmO to be a rather unique system. The similarities of the 
crystal symmetry and the lattice constant of SmO to the well studied ferromagnetic semiconductor 
EuO, makes SmO/EuO thin film interfaces an excellent contender towards realizing the quantum 
anomalous Hall effect in a strongly correlated electron system. 


Topological insulators (TI) are materials which exhibit 
a fundamentally new physical phenomena that was first 
predicted by theorists [mamzi and subsequently ascer¬ 
tained in experiments Although most work con¬ 

siders band semiconductors, the concept of topology can 
be extended to correlated systems for which many excit¬ 
ing new effects can be expected. In analogy to correlation 
driven fractional quantum Hall states one might expect 
fractional Chern insulating states to emerge in correlated 
topological insulators [T3HI7]. The experimental realiza¬ 
tion of such states would open a whole new field of pos¬ 
sibilities. 

SmBe was recently reported to be a correlated mixed 
valent topological insulator [TBHIS] . The highly corre¬ 
lated flat Sm-/ derived bands hybridize with the disper¬ 
sive and itinerant Sm-d derived bands to form a mixed 
valent insulating ground state. Although the insulat¬ 
ing state and topology of the bulk of SmBe is well de¬ 
fined, [001] surface of SmBe is polar and therefore in¬ 
herently unstable [24] . This hinders the unique determi¬ 
nation of topological surface states and in turn inhibits 
the creation of robust technological devices. Here, us¬ 
ing density functional theory (DFT) based calculations, 
we show that SmO is a mixed valent correlated com¬ 
pound with a band structure similar to the topological 
nontrivial high-pressure gold phase of SmS [25|. In con¬ 
trast to SmS, SmO is predicted to have a topological 
nontrivial ground state at ambient pressure. Addition¬ 
ally, our calculations show that the topological nontrivial 
ground-state of SmO is stable for a large range of Sm- 
O distances, including both positive and negative strain. 
The Sm-/ band filling is thus tunable by strain, which 
opens up the possibility to create correlated topological 
nontrivial bands at different filling. Subsequent to the 
enormous success in the design and fabrication of semi¬ 
conductors, it is presently possible to grow high-quality 
oxide thin films and heterostructures PHUT]. SmO thus 
seems the ideal candidate for the experimental realiza¬ 


tion of topological nontrivial correlated flat bands. In ad¬ 
dition, experimental realization of the theoretically pre¬ 
dicted quantum anomalous Hall effect (QAHE) (quan¬ 
tized Hall conductance in the absence of an external 
magnetic field) has been challenging due to difficulties 
in obtaining insulating bulk TTs concurrent with homo¬ 
geneous magnetic doping m Recently, notwithstanding 
the presence of bulk carrier density (i.e. not a nominally 
insulating bulk), long range ferromagnetism and QAHE 
was successfully observed for the first time in Cr-doped 
(Bi,Sb) 2 Te 3 TI thin films [30] . Other proposals to¬ 
wards realizing a QAHE are to either grow a TI on top 
of a ferromagnetic insulator iniiss; transition-metal ox¬ 
ide heterostructures, [33l|34] or to deposit a layer of heavy 
atoms with large spin-orbit coupling on a magnetic insu¬ 
lator [35] . Fortuitously, a well studied substrate with the 
same symmetry as SmO is readily available in the form 
of the ferromagnetic semiconductor EuO [36ll42] . Con¬ 
sequently, we propose growing SmO/EuO thin film in¬ 
terfaces to test the feasibility of obtaining a QAHE in a 
strongly correlated electron system. 

SmO crystallizes in the rock-salt structure with 
a = 4.9414-4.943 A ga ig. Considering the lattice pa¬ 
rameters of trivalent neodymium chalcogenides and inter¬ 
polating between neodymium and terbium compounds, 
Leger et. al. obtain a lattice constant of 4.917 A for 
Sm^+0[43]. Similarly, considering other divalent lan¬ 
thanide oxides, a lattice constant of 5.15 A is expected 
for Sm^+0. Ergo, the experimental cell parameter of fee 
SmO lies between that of Sm^+0 and Sm^+0, which lead 
the authors [43] to conclude that samarium is in an in¬ 
termediate valence state in SmO. Using Vegard’s law, 
a valence of 2.92 was assigned [43] . Electrical resistiv¬ 
ity measurements as a function of temperature reveals 
a like increase between 4.2 and 20 K, followed by a 
rapid increase up to 32 K. The resistivity shows a linear 
temperature dependence above 60 K [45] . Krill and co¬ 
workers [44] measured the magnetic susceptibility of SmO 
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FIG. 1. (Color online) Phase diagram (d) of SmO as a function of the fee lattice parameter within LDA+SO+f/ {U = 6 eV, 
Jh = 0 eV) along with representative FPLO band structures for the (a) nontrivial d — f band inversion, (b) trivial insulator and 
(a) nontrivial s — f band inversion scenarios. The 4/, 5d and 6s orbital character derived bands are represented by red, blue 
and green colored symbols respectively. The size of the symbols represent the weight of the various orbital contributions to 
the underlying bands. The Brillouin zone of an fee lattice is displayed along with the eight time reversal invariant momenta 
(TRIM) points. 


and observed the roughly constant magnetic susceptibil¬ 
ity and the non-divergence of the low-temperature sus¬ 
ceptibility. They drew parallels to the susceptibilities of 
other intermediate valence compounds SmBe and gold- 
SmS under pressure, wherein the low temperature behav¬ 
ior is not explainable by crystal field effects alone [44]. 

The electronic structure calculations are performed us¬ 
ing the full-potential non-orthogonal local orbital code 
(FPLO) [46]. The local density approximation (LDA) with 
the Perdew and Wang flavor m of the exchange and cor¬ 
relation potential was chosen. To account for the strong 
spin-orbit (SO) coupling of the / electrons, we employ 
the full four-component relativistic scheme. Additionally, 
the strong Coulomb repulsion between the 4/ electrons of 
samarium are included in a mean-field way by applying 
LDA+SO+[/ with the ’’fully localized limit” (FLL) dou¬ 
ble counting term [48] . This level of theory will not suffice 
to describe quantitatively the Sm 4/ spectral weight and 
its dispersion, but is sufficient to predict the Sm d — f 
inversion. We have tested the robustness of these pre¬ 


dictions by varying U (5 to 7 eV) and Jr (0 to 0.7 eV) 
and by using different functionals (LDA and GGA us¬ 
ing FPLO and the modified Becke-Johnson approach [49] 
with [7=3 eV and Jh=0 eV using Wien2k[50]). Impor¬ 
tantly, as pointed out by Martin and Allen [2] for SmBe 
and related systems, the symmetry of the relevant Sm 
one particle orbitals and the many body electron removal 
Green’s function is the same. This allows for a possible 
adiabatic continuation from the LDA results as presented 
here to the full interacting system without changing the 
topology of the system. See Supplemental Material [52] 
for further discussion on the Sm 4/ many body state. 

Collected in Fig.|Sl[a) is the FPLO non-spin polarized, 
full-relativistic band structure of SmO with the inclusion 
of the strong Coulomb interaction (LDA+SO+7/) for the 
experimental lattice constant. With samarium being in 
the 2+ configuration, the 4/ states are split into lower 
lying and filled 4 / 5/2 states that can accommodate 6 elec¬ 
trons and higher lying (above 5.5 eV), empty 4 / 7/2 states 
(see Supplemental Material). We do not observe a direct 
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FIG. 2. (Color online) Left: Surface Bloch spectral density 
of the first 12 SmO layers of a semi-infinite solid 
with [001]-surface. Note, that there are two Dirac cones at 
the X-point, while the one at the F-point falls into the bulk 
projected band structure and hence forms a surface resonance 
(see right panel). Right: AB\{k) with downwards shifted 4/- 
electron pocket around the f-point to reveal the third Dirac 
cone. 


energy gap around the Fermi level (Ep). The material is 
semi-metallic with a “warped gap” such that everywhere 
in the BZ (at Ep), band number x is always below band 
number x+1. There are small electron and hole pockets 
at X and F respectively, but there are no band crossings 
between the highest occupied 4 / 5/2 bands and the low¬ 
est unoccupied bands. Besides, varying U (5 to 7 eV) or 
Jh (0 to 0.7 eV) does not change the above mentioned 
features, since the oxygen 2 j 9 , samarium 5d and 65 bands 
experience a constant shift with respect to the localized 
4/ states. There exists only one report m in literature 
on the transport properties of SmO wherein the resistiv¬ 
ity decreases with decreasing temperature and thus could 
hint to metallic character, consistent with the lack of an 
energy gap in our calculations. In contrast to the local¬ 
ized and not very dispersive 4/ bands, the 5d bands of 
samarium are dispersive, broad and dip below the Fermi 
level, retaining a 100% weight at the X point. As a conse¬ 
quence, we observe a 4/-5d band inversion at the X point, 
resulting in a nontrivial topology, similar in sense to the 
band inversions reported for SmBe and the high pressure 
gold phase of SmS. This band inversion and the resulting 
topological indices 1;(000) allows us to classify SmO as 
a 3D strongly topological semi-metal (more details are 
provided in Supplemental Material). 

A well known issue with LDA when dealing with semi¬ 
conductors, is the underestimation of band gaps. Since 
nontrivial topology depends on band inversions, under¬ 
estimation of band gaps could sometimes lead to wrong 
sequence of the underlying band structure. For example, 
in HgTe, a s-p TI, LDA correctly predicts the band in¬ 
version, but results in an incorrect band sequence at the 


FIG. 3. (Golor online) Maximum possible / valence due to 
/ —d 01 f —s promotion in SmO while retaining the 
nontrivial topology as a functi on of lattice parameter. The 
color map follows that of Fig.|Sl[d): light blue = nontriv¬ 
ial topology from d — f inversion, light green = nontrivial 
topology from s — f inversion, light pink = trivial topology. 
The equilibrium lattice parameters for a high-spin 3+ and 2+ 
state are computed using L(S)DA+SO+f/ scheme (dashed red 
line). The estimates from an empirical extrapolation for the 
limiting cases are obtained from Ref. 23] (dashed green line). 
The filled green circle denotes the experimental bulk SmO. 
SmO remains a topological semi-metal for a wide range of ten¬ 
sile strain (> 1%), since the calculated maximum valence, is 
above both the empirical extrapolation and L(S)DA+U+SO 
estimates. 


F point. The correct band sequence is obtained by using 
the MBJLDA exchange potential [53]. Literature on the 
efficiency of MBJLDA over the traditional LDA+U for 
/ electron systems is still scarce. In a recent work on 
SmS, the authors have employed MBJLDA+SO+t/ {U 
= 3 eV) to open a band gap of 0.2 eV in the ambient 
pressure black phase such that the gap is consistent with 
the available experiments |25|. Using the same param¬ 
eters, the authors conclude that the high pressure gold 
phase of SmS is a topological metal. We compared the 
LDA+SO+U results with that of MBJLDA+SO and ob¬ 
tain a consistent picture for the band inversion between 
the two approaches (see Supplemental Material). No di¬ 
rect band gap is opened with MBJLDA+SO+U (U = 3 
eV), in accordance with our LDA+SO+U results and as 
well as with the experimental report m- 

To provide additional confirmation of the topological 
non-triviality in SmO and to explicitly identify the pro¬ 
tected surface states, we have calculated the Bloch spec¬ 
tral density {A-B\{k)) of the 12 topmost surface layers of 
a semi infinite solid [54| with [001] surface. We obtain 
two weakly interacting quasi-degenerate Dirac-cones at 
M, and a single Dirac-cone at F which is hidden in the 
bulk and becomes a surface resonance due to the semi- 
metallic nature of SmO. To clearly identify the Dirac- 
cone at F, we have repeated the slab calculations such 
that the hole pockets at F for the bulk are pushed down. 
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This in turn, allows one to unambiguously identify all 
three Dirac-cones in a [001] terminated SmO. More de¬ 
tails are provided in the Supplemental Material. 

Having established the nontrivial topology in SmO for 
the experimental bulk lattice parameter, we consider the 
scenario of growing thin films of SmO. In general, the 
lattice parameter of the substrate plays a decisive role in 
determining the lattice parameter of the thin film. Then, 
the relevant question to answer is the robustness of the 
topological semi-metal state as a function of lattice pa¬ 
rameter variation. To this end, we have investigated the 
topological indices for various lattice parameters. We 
have retained the cubic symmetry of the unit cell (i.e 
isotropic volume change) in SmO, based on the experi¬ 
mental studies on isostructural samarium systems, which 
evidences a negative value of the Poisson ratio v (elastic 
constant Cyi < 0 ), characteristic of an isotropic volume 
change [55l |56]. Since the nontrivial topology is obtained 
due to d — f band inversion at the X point, SmO is topo¬ 
logical as long as the 5d band bottom is below the Ep. 
The nontrivial topology is suppressed when the 5d band 
bottom moves above Ep or when other bands dip below 
Ep (Eig.jS^b)). In the case of SmO, it is the 6s band 
bottom at the T point that shifts to lower values as a 
function of increasing lattice parameter. The shift of the 
5d and 6s band bottoms (at X and T respectively) with 
respect to the Ep as a function of lattice parameter are 
collected in Eig.j^^d) for LDA+SO+t/ approach. Treat¬ 
ing the strong 4/ correlations on a mean-field level, the 
topological semi-metal state is quite robust and remains 
so, up to 5.38A, a 9% increase in lattice parameter. On 
an experimental level, this result is very promising, since 
it provides a large range of substrate lattice constants for 
which SmO thin films reveal topological non-triviality. 
During the evaluation of the topological indices using 
LDA+SO+I/, we became aware of another interesting 
feature, a nontrivial topology due to s — f band inversion 
at T for 11% and larger lattice constants (> 5.5 A). In 
LDA+SO+I/, the 6d band bottom shifts above the Ep 
and becomes trivial before the 65 band bottom dips be¬ 
low Ep. Eurther increase in the lattice parameter then 
results in a 5 — / band inversion at T which is again 
topological (Eig.|SI[c)). Although, the probability for a 
successful growth of SmO on a substrate with a 11% and 
larger lattice constant may be low, it nonetheless offers 
another exciting route towards realizing a topologically 
nontrivial state, only this time with as — / band inver¬ 
sion. We verified the validity of our results for a range 
of U values ( 5 to 7 eV) and as well as for the MBJLDA 
exchange potential. 

In the calculations so far, we have allowed Sm to be in 
the intermediate valence state and have treated this on 
a DET level. The mixed valent situation can be viewed 
as a 4/ to 6d promotion of a certain fractional amount 
of electrons due to the hybridization of the localized 4/ 
states with the band like 5d states. Based on the handful 


of experimental reports available in literature, the va¬ 
lence assigned to Sm in SmO is 2.92, i.e. a promotion 
of 0.92 electrons. To actually calculate the amount of 
electron promotion is a difficult task, since DET in the 
Kohn-Sham scheme is based on a single Slater determi¬ 
nant approximation and prevents a correct description of 
a many body intermediate valence state. In other words, 
the position of the 4/ bands relative to the 6d/6s as 
calculated in DET is not a reliable quantity (see Sup¬ 
plemental Material). Yet, we can address the effect of 
mixed valency on the topological properties of SmO in 
a rigorous quantitative way, by calculating the amount 
of electrons (integration of the density of states) that 
can be accommodated in the dispersive 6d or 6s band, 
thereby mimicking various amounts of 4/ — 6d or 4/ 
— 6s promotion. Eor the lattice constants which has the 
band inversion at the X point, we calculate the amount 
of electrons that can be contained in the 5d band be¬ 
fore beginning to populate the 6s band. Note that the 
nontrivial topology is maintained as long as only the 5d 
is occupied and the 6s remains unoccupied. Eig.[^ dis¬ 
plays this electron amount as a function of the lattice 
constant. Analogously, on the other side of the phase 
diagram, wherein non-triviality is manifested due to a 
s — f inversion, we calculate the amount of electrons 
that can be contained in the 6s band before beginning 
to populate the 5d band. Equilibrium lattice constants 
for the 3+ and 2+ limiting cases are estimated using 
spin-polarized L(S)DA+SO+I/ schemes (see Supplemen¬ 
tal Material). We observe that the theoretical estimated 
valence and lattice constants always yield a topological 
nontrivial ground state. Though the estimate from an 
empirical extrapolation for the SmO valence places the 
bulk in the trivial region of the phase diagram, nontrivial 
topology is quickly reinstated for lattice parameters with 
a small tensile strain (> 1%). Spectroscopy experiments 
are necessary to confirm the samarium valency, and posi¬ 
tions of the samarium s and d bands with respect to the 
/ states in bulk SmO and SmO thin films. 

By virtue of the large surface to volume ratio, thin films 
of systems with nontrivial topology with dominant sur¬ 
face states are highly sought-after. Owing to the simple 
fee symmetry and the flexibility in observing the topo¬ 
logical ground state in SmO for a wide range of lattice 
parameters, we anticipate plenty of options for suitable 
substrates. One particular substrate that invokes special 
interest is the fee ferromagnetic semiconductor EuO 
1^ . which has a lattice constant of 5.14 A, only a 4% lat¬ 
tice mismatch with that of bulk SmO. EuO by itself is 
an attractive material with many functionalities, includ¬ 
ing metal-insulator transition, magnetic phase transition, 
colossal magneto-resistance, etc[39H42]. These function¬ 
alities have produced an abundance of research on grow¬ 
ing high quality thin films of EuO. So, a quick progress 
in the engineering of SmO/EuO thin film interfaces can 
be expected. Additionally, this would open up the possi- 
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bility for another experimental realization of the QAHE, 
this time in a strongly correlated electron system with 
topologically nontrivial flat bands. Despite the lack of a 
bulk band gap in Cr-doped (Bi,Sb) 2 Te 3 thin films, QAHE 
was observed at 30 mK Enmu]. In a similar manner, the 
presence of a semi-metallic bulk in SmO will give rise to 
bulk carriers, but given the tiny electron and hole pockets 
around X and E respectively, and by making use of thin 
films, the amount of bulk charge carriers can be kept 
quite small and should not completely overshadow the 
QAHE. 

In summary, we have investigated in detail the elec¬ 
tronic structure of SmO and established the nontrivial 
topology of the band structure such that SmO can be 
classified as a 3D strongly topological semi-metal. The 
topological state in SmO is very robust and prevails for 
a wide range of lattice parameters, which leads us to 
propose SmO as an ideal candidate for further investiga¬ 
tions of topologically nontrivial correlated flat bands in 
bulk and thin-film form. 

The authors thank Prof. J. Allen, S. Wirth and P. 
Thalmeier for discussions. D.K. acknowledges funding 
by the DEG within EOR 1346. 


[1] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 
(2005). 

[2] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 
(2005). 

[I] L. Fu and C. L. Kane, Phys. Rev. B 76, 045302 (2007). 

[4] J. E. Moore and L. Balents, Phys. Rev. B 75, 121306(R) 
(2007). 

[5] R. Roy, Phys. Rev. B 79, 195321 (2009). 

[6] L. Fu, C. L. Kane, and E. J. Mele, Phys. Rev. Lett. 98, 
106803 (2007). 

[7] H. Zhang, C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang, and S.-C. 
Zhang, Nat. Phys. 5, 438 (2009). 

[8] M. Konig, S. Wiedmann, C. Brune, A. Roth, H. Buhmann, 
L. W. Molenkamp, X. L. Qi, and S. C. Zhang, Science 
318, 766 (2007). 

[9] D. Hsieh, D. Qian, L. Wray, Y. Xia, Y. S. Hor, R. J. Cava, 
and M. Z. Hasan, Nature 452, 970 (2008). 

[10] D. Hsieh, Y. Xia, L. Wray, D. Qian, A. Pal, H. Dil, J. Os- 
terwalder, F. Meier, G. Bihlmayer, C. L. Kane, Y. S. Hor, 

R. J. Cava, and H. M. Z., Science 323, 919 (2009). 

[II] Y. Xia, D. Qian, D. Hsieh, L. Wray, A. Pal, H. Lin, 
A. Bansil, D. Grauer, Y. S. Hor, R. J. Cava, and M. Z. 
Hasan, Nat. Phys. 5, 398 (2009). 

[12] Y. L. Chen, J. G. Analytis, J.-H. Chu, Z. K. Liu, S.-K. 
Mo, X. L. Qi, H. J. Zhang, D. H. Lu, X. Dai, Z. Fang, 

S. C. Zhang, 1. R. Fisher, Z. Hussain, and Z.-X. Shen, 
Science 325, 325 (2009). 

[13] E. Bergholtz and Z. Liu, Int. J. Mod. Phys. B 27, 1330017 
(2013). 

[14] Z. Liu, E. J. Bergholtz, and E. Kapit, Phys. Rev. B 88, 
205101 (2013). 

[15] S. Flach, D. Leykam, J. D. Bodyfelt, P. Matthies, and 
A. S. Desyatnikov, Euro Phys. Lett. 105, 30001 (2014). 


[16] S. Kourtis, T. Neupert, C. Chamon, and C. Mudry, Phys. 
Rev. Lett. 112, 126806 (2014). 

[17] A. G. Grushin, A. Gomez-Leon, and T. Neupert, Phys. 
Rev. Lett. 112, 156801 (2014). 

[18] M. Dzero, K. Sun, V. Galitski, and P. Coleman, Phys. 
Rev. Lett. 104 , 106408 (2010). 

[19] M. Dzero, K. Sun, P. Coleman, and G. V., Phys. Rev. 
B 85, 045130 (2012). 

[20] T. Takimoto, J. Phys. Soc. Jpn. 80, 123710 (2011). 

[21] S. Wolgast, t. Kurdak, K. Sun, J. W. Allen, D.-J. Kim, 
and Z. Fisk, Phys. Rev. B 88, 180405(R) (2013). 

[22] D. J. Kim, J. Xia, and Z. Fisk, Nature Mat. 13, 466 
(2014). 

[23] S. Suga, K. Sakamoto, T. Okuda, K. Miyamoto, 
K. Kuroda, A. Sekiyama, J. Yamagichi, H. Fujiwara, 
A. Irizawa, T. Ito, S. Kimura, T. Balashov, W. Wulfhekel, 
S. Yeo, F. Iga, and S. Imada, J. Phys. Soc. Jpn. 83, 014705 
(2014). 

[24] Z.-H. Zhu, A. Nicolaou, G. Levy, N. P. Butch, P. Syers, 
X. F. Wang, J. Paglione, G. A. Sawatzky, 1. S. Elfimov, 
and A. Damascelli, Phys. Rev. Lett. Ill, 216402 (2013). 

[25] Z. Li, J. Li, P. Blaha, and N. Kioussis, Phys. Rev. B 89, 
121117(R) (2014). 

[26] J. Herber, Nature 459, 28 (2009). 

[27] J. Mannhart and D. G. Schlom, Science 327 , 1607 (2010). 

[28] H. Ji, J. M. Allred, N. Ni, J. Tao, M. Neupane, A. Wray, 

S. Xu, M. Z. Hasan, and R. J. Cava, Phys. Rev. B 85, 
165313 (2012). 

[29] C. Chang and et. al, Science 340 , 167 (2013). 

[30] H. Ke, M. Xu-Cun, C. Xi, L. Li, W. Ya-Yu, and X. Qi- 
Kun, Chin. Phys. B 22, 067305 (2013). 

[31] L. Fu and C. L. Kane, Phys. Rev. Lett. 102, 216403 
(2009). 

[32] A. R. Akhmerov, J. Nilsson, and C. W. J. Beenakker, 
Phys. Rev. Lett. 102, 216404 (2009). 

[33] A. M. Cook and A. Paramekanti, Phys. Rev. Lett. 113 , 
077203 (2014). 

[34] T. Cai and et at, arXiv:1310.2471 (2013). 

[35 K. F. Garrity and D. Vanderbilt, Phys. Rev. Lett. 110 , 
116802 (2013). 

[36] B. T. Matthias, R. M. Bozorth, and J. H. Van Vleck, 
Phys. Rev. Lett. 7 , 160 (1961). 

[37] W. Nolting, Phys. Stat. Sol. (B) 96, 11 (1979). 

[38] G. Adachi and N. Imanaka, Chem. Rev. 98, 1479 (1998). 

[39] P. G. Steeneken, L. H. Tjeng, 1. Elfimov, G. A. Sawatzky, 
G. Ghiringhelli, N. B. Brookes, and D. J. Huang, Phys. 
Rev. Lett. 88, 047201 (2002). 

[40] A. Schmel, V. Vaithyanathan, A. Herrnberger, S. Thiel, 
C. Richter, M. Liberati, T. Heeg, M. Rockerath, L. Kourk- 
outis, S. Miihlbauer, P. Boni, D. A. Muller, Y. Barash, 
J. Schubert, Y. Idzerda, J. Mannhart, and D. Schlom, 
Nat. Mater. 6, 882 (2007). 

[41] R. Sutarto, S. G. Altendorf, B. Coloru, M. Moretti Sala, 

T. Haupricht, C. F. Chang, Z. Hu, C. Schiissler- 
Langeheine, N. Hollmann, H. Kierspel, H. H. Hsieh, L. H. 
J., C. T. Chen, and L. H. Tjeng, Phys. Rev. B 79, 205318 
(2009). 

[42] A. Melville, T. Mairoser, A. Schmehl, M. Fischer, 
S. Gsell, M. Schreck, D. D. Awschalom, T. Heeg, B. Hol¬ 
lander, J. Schubert, and D. G. Schlom, App. Phys. Lett. 
103, 222402 (2013). 

[43] J. M. Leger, N. Yacoubi, and J. Loriers, Inorg. Chem. 
19 , 2252 (1980). 

[44] G. Krill, M. F. Ravet, J. P. Kappler, L. Abadli, J. M. 



6 


Leger, N. Yacoubi, and C. Loriers, Solid State Comm. 
33, 351 (1980). 

[45] J. M. Leger, P. Almonino, J. Loriers, P. Dordor, and 
B. Coqblin, Phys. Lett. 80A, 325 (1980). 

[46] K. Koepernik and H. Eschrig, Phys. Rev. B 59, 1743 
(1999). 

[47] J. P. Perdew and Y. Wang, Phys. Rev. B 45, 13244 
(1992). 

[48] M. T. Czyzyk and G. A. Swatzky, Phys. Rev. B 49, 14211 
(1994). 

[49] F. Tran and P. Blaha, Phys. Rev. Lett. 102 , 226401 
(2009). 


[50] P. Blaha and K. Schwarz, Comput. Mater. Sci. 28, 259 
(2003). 

[2] R. Martin and J. Allen, J. Appl. Phys. 50, 7561 (1979). 

[52] See Supplemental Material at 

http://link.aps.org/Supplemental/ .. 

[53] H. Zhang and S. Zhang, Phys. Status Solid! RRL 7, 72 
(2013). 

[54] J. Velev and W. Butler, J. Phys.: Condens. Matter 16, 
R637 (2007). 

[55] B. Liithi, J. Magn. Magn. Mat. 52, 70 (1985). 

[56] P. Wachter, Handbook on the Physics and Chemistry of 
Rare Earths (Vol. 19, Chapter 132, 1998). 


Supplemental Material 


CALCULATION OF TOPOLOGICAL INDICES 

To calculate the Z 2 topological index for a system with 
inversion symmetry, we use the parity criteria as pro¬ 
posed by Fu and Kane[T|, wherein the product of the 
parities of all occupied Kramers doublets at each time re¬ 
versal invariant momentum (TRIM) is determined. For 
a fee lattice, and hence a bee Brillouin zone, we have the 
following 8 TRIM points : 1 T (0,0,0); 3 X (l,0,0)27r/a 



FIG. SI. (Color online) The calculated FPLO band struc¬ 
ture of SmO (a = 4.941 A) within a non-spin polarized, full- 
relativistic approximation including strong Coulomb correla¬ 
tion (LDA+SO+U). A U and Jh value of 6 eV and 0 eV have 
been used respectively. The size of the symbols represent the 
weight of the various orbital contributions to the underlying 
bands. The numbers within parentheses are band indices to 
assist in assigning the various topological indices in Table. [ST] 


and 4 L (|,|,|)27r/a. The sign of the product of the pari¬ 
ties of all occupied doublets (Kramers degenerate bands) 
at these eight TRIM points refer to a trivial topology, 
when positive and to a nontrivial topology, when neg¬ 
ative. Under spatial inversion, the d orbitals are even 
while the / orbitals are odd. From the calculated band 
structures, it is clear that the occupied bands between 
0 to -2 eV at F and the four L points possess / orbital 
character, while due to the 4/-5d band inversion, one of 
the occupied band at the three X points possesses d or¬ 
bital character. The product of the parities for all the 
occupied bands will then be negative, giving rise to a 
nontrivial topology in the band structure of SmO. This 
result is valid for both LDA+SO and LDA+SO+U calcu¬ 
lations. This qualitative discussion was confirmed by ex¬ 
plicitly calculating the space group representations of the 
inversion operator of all bands at the TRIM points (Ta- 
ble|^. This gives access to the four topological indices 
^ 0 ; (^i^ 2 ^ 3 )- The resulting indices for all bands up to the 
highest occupied 4 / 5/2 bands are 1; (000), which makes 
SmO strongly topological. As a result of the particular 
arrangement of the bands around the Ep in SmO, which 
creates a semi-metallic “warped gap” (no band crossings 
at Ep), it is possible to apply the parity counting scheme 
at the TRIM points, since this “warped gap” ensures the 
continuity condition for the phases of the wave functions 
in the whole BZ which leads to the counting scheme. 
In consequence one could classify SmO as a 3D strongly 
topological semi-metal. 

The topological indices switch from trivial to nontrivial 
at the bottom of the six 4 / 5/2 derived bands and switch 
back to trivial above the band, which forms the elec¬ 
tron pockets around the X point, which in consequence 
tells us that topological surface bands could bridge en¬ 
ergy window from the lowest 4 / 5/2 band to the lowest 
unoccupied 5d band. In Table we provide the details 
of the topological character of the important bands. The 
energy separation of about 2 eV between the 02p bands 
and the lowest 4/ state ensures that the gap between 
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TABLE SI. Calculated topological indices of selected bands 
in SmO within LDA+SO+f/ approximation (a = 4.941 A). 
Sn refers of the band energy at T; n refers to the band in¬ 
dices, shown in Fig. |Sl| T, X and L are the TRIM points. 
vq] {iyiU 2 iy 3 ) are the four topological indices. The parity of 
the Kramers doublet at the various TRIM points are denoted 
by O (odd) and E (even). Additionally, the dominating or¬ 
bital character of the bands at the TRIM points are listed. 


Sn (r) eV, (n) 

r 

3xX 

4xL 

z/o; (1^11^21^3) 

-3.650 

0 (2p) 

0 (2p) 

E (2p) 

0; (000) 

-0.493 (1,2) 

0 ( 4 / 5 / 2 ) 

E (5d) 

0 ( 4 / 5 / 2 ) 

1;(000) 

+0.030 (3,4) 

0 ( 4 / 5 / 2 ) 

0 ( 4 / 5 / 2 ) 

0 ( 4 / 5 / 2 ) 

1; (000) 

+0.030 (5,6) 

0 ( 4 / 5 / 2 ) 

0 ( 4 / 5 / 2 ) 

0 ( 4 / 5 / 2 ) 

1; (000) 

+1.457 (7,8) 

E (6s) 

0 (5d) 

E (5d) 

0; (000) 


02p and Sm 4/ must be trivial due to the fact that SO 
coupling is too small to create band inversions bridging 
this gap. Hence the highest 02p band forms the trivial 
base line for discussing the topological properties of the 
remaining bands. In the table, we list the energy of the 
bands at the T point, the band indices (as indicated in 
Fig.[ST]), the parity of the Kramers doublets (O: odd, E: 
even) together with the dominating orbital character at 
the three symmetry distinguished classes of TRIM points 
as well as the resulting topological indices obtained from 
all bands including the ones mentioned in the first col¬ 
umn. Note that, bands numbered 3.. .6 are trivial bands, 
and do not change the topology. 


BAND STRUCTURE USING MODIFIED BECKE 
JOHNSON APPROACH 


To make sure that the topological characteristics are 
not dependent on the choice of exchange and correlation 
functional, we calculated the band structure of SmO us¬ 
ing the recently proposed modified Becke-Johnson func¬ 
tional (MBJLDA) as implemented in WIEN2K. Fig. S2 
shows the band structure obtained using MBJLDA+SO. 
The f — d band inversion at X point is consistent to our 
LDA+SO+U calculations, though the hybridization gap 
and the warped band-gap in MBJLDA+SO are smaller 
(but, of the same order). 


SURFACE STATE CALCULATIONS 

To asses the nature of the surface band structure we 
calculated the Bloch spectral density {AB\{k)) of the 12 
topmost surface layers of a semi infinite solid with [001] 
surface. For this purpose the DFT band structure was 
fitted with atom centered maximally projected Wannier 
functions of Sm Sd, 6s and 4/ and O 2p character. A 
soft energy cutoff was employed to achieve band disen¬ 
tanglement at +6eV. The resulting fit models all relevant 




FIG. S2. (Color online) Band structure of SmO with or¬ 
bital character using MBJLDA+SO approximation as imple¬ 
mented in WIEN2K. The topology and order of the bands are 
consistent to the ones obtained using LDA+SO+U scheme. 
The blow-up in the bottom panel clearly shows the opening of 
the warped band-gap along T-X, but smaller in size compared 
to FPLO. For clarity we have plotted the bands without any 
band character. 


bands around the Fermi level within a few meV accuracy. 
Afterwards, the obtained bulk hopping parameters were 
mapped onto a semi infinite solid and a Greens function 
technique was used to calculate AB\{k) for the surface lay¬ 
ers. The resulting surface bands clearly show two weakly 
interacting Dirac cones around the surface projected M- 
point. In order to be of topology induced nature, there 
must be an odd number of Dirac cones. The surface Bril- 
louin zone for a [001] surface ends up with two bulk X- 
points being projected onto the surface M-point and one 
bulk X-point being projected onto the surface f-point 
(see Fig. S3). Hence, one expects two Dirac cones at M 
and one at F. The back-folding of the bulk band struc¬ 
ture due to projection onto [001] leads to the appearance 
of electron pocket states from the bulk X-point at f, 
which overlap with the bulk projected states of the bulk 
r hole pocket at T. Hence, the third Dirac cone will be 
immersed into projected bulk states and can only form a 
surface resonance. In order to prove the existence of this 
third Dirac cone, we performed an additional calculation. 





































where the 4/ electron pocket around the bulk F-point 
are lowered in energy by application of a k-dependent 
potential. This removes bulk projected states from the 
low energy region at the surface F point and reveals the 
third Dirac cone. Finally, to rule out drastic effects of 
the relaxation of the surface electronic structure we per¬ 
formed a full DFT slab calculation for 21 SmO layers. 
The resulting slab band structure shows the same two 
weakly interacting Dirac cones around the M-point as 
the semi-infinite calculation. Slight shifts of bands oc¬ 
cur, but the overall structure equals that of the mapped 
model calculation. 

NOTES ON THE MANY BODY STATE OF SMO 

The Sm 4/ states are known to be strongly correlated. 
Nonetheless one can obtain useful information about this 
system using DFT within the local density approxima¬ 
tion. The large onsite Coulomb repulsion between the / 
electrons restricts the local valence occupations to fluc¬ 
tuate between the and configurations. The lattice 
constants for the limiting cases of a pure 2+ (/^) and 
pure 3+ (/^) SmO compound are determined using spin 
polarized L(S)DA+SO+U. The local atomic (/^) con¬ 
figuration has a lowest Hunds-rule multiplet state which 
belongs to the ^F term. LDA in the Kohn-Sham 

scheme does capture these local states as they are sin¬ 
gle Slater determinant representable (as any Hunds-rule 
high spin ground-state). 

In order to determine the topology of the bands and 
in order to answer the question if hybridization between 
the Sm d derived bands and the Sm / derived bands is 
allowed, it is important to have the local symmetry cor¬ 


rect. The spin polarized state found in L(S)DA does not 
represent the local Sm character correctly. Within 
L(S)DA the Sm configuration has a local moment 
{Mz = 2Sz + I/Z = 3) whereas the many body ground- 
state is a singlet ((J^) = J{J + 1) = 0) without a local 
moment [Lz = 0^ Sz = 0, Mz = 0). The many-body 
ground-state without a local moment is consistent with 
experimental observations on divalent Sm compounds. In 
order to reproduce the local symmetry of the many body 
state of an configuration (^Fq term with J = 0, be¬ 
longing to the Aig representation) we use non-spin polar¬ 
ized DFT calculations. In this case the configuration 
is represented by a state with six electrons in the j = 5/2 
bands. Although this state has an incorrect expectation 
value of L ((T^) = 24/7) and S {{S‘^) = 24/7) com¬ 
pared to the many body state (and thus a to high local 
Coulomb energy), it does have the right symmetry (^i^) 
and total momentum (J = 0). The configuration is 
represented by single hole excitations starting from the 
configuration. A peculiarity of starting from a many 
body state with Aig symmetry is that the many body 
one electron excitations have the same symmetry as the 
independent electron bands found in LDA, as shown by 
Martin and Allen[2]. 
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FIG. S3. (Color online) Bulk and surface Brillouin zone of a 
fee lattice. 






